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Abstract 



The rules of soft-collinear effective theory can be used naively to write hard scattering cross-sections as convolutions of separate 
hard, jet, and soft functions. One condition required to guarantee the validity of such a factorization is the infrared safety of these 
functions in perturbation theory. Using e + e~ angularity distributions as an example, we propose and illustrate an intuitive method to 
<3\ test this infrared safety at one loop. We look for regions of integration in the sum of Feynman diagrams contributing to the jet and 
soft functions where the integrals become infrared divergent. Our analysis is independent of an explicit infrared regulator, clarifies 
how to distinguish infrared and ultraviolet singularities in pure dimensional regularization, and demonstrates the necessity of taking 
zero-bins into account to obtain infrared-safe jet functions. 
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1. Introduction 

C"| Factorization restores predictive power to calculations in 
Q_i Quantum Chromodynamics (QCD) which cannot be carried out 
Q^exactly due to the contributions of nonperturbative effects. By 
separating perturbatively-calculable and nonperturbative con- 
tributions to observables in QCD and relating the nonpertur- 
bative contributions to different observables to each other, we 
^vq gain the ability to make real predictions. 
J> Proving factorization rigorously is a technically challeng- 
I" ing undertaking, which traditionally has been formulated in full 

QCD din. More recently, many formal elements of these fac- _ 1 yi ^ 

I torization proofs, such as power counting, gauge invariance, the a Q ^ ' 
• organization of soft gluons into eikonal Wilson lines, and their 
decoupling from collinear modes, have been organized in the 
framework of soft-collinear effective theory (SCET) 



An ideal set of observables for which to examine factoriz- 
ability is the set of angularities r a lfl3l . which are two-jet e + e~ 
event shapes dependent on a tunable parameter a controlling 
how sensitive the event shape is to radiation along the jet axes 
or at wider angles. Varying a between and 1 interpolates be- 
tween the thrust lfT4l[T31 and jet broadening 1 16 1 event shapes, 
but a can take any value — oo < a < 2 and give an infrared-safe 
observable in QCD. Angularities are known to be factorizable, 



however, only for a < 1 [13|. For events e + e 
larity of a final state X is 



X, the angu- 
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These generic properties of the effective theory allow one to 
write at least nominally a formula "factorized" into collinear 
(jet) and soft functions for an arbitrary hard scattering cross- 
section in which strongly-interacting light-like particles partic- 
^ ipate Q. Examples are the factorization of a large class of two- 
jet event shape distributions in e + e~ annihilations to light quark 
jets |8, 9 10 1, jet mass distributions for e + e~ to top quark jets 
ifTTIl . or arbitrary jet cross-sections in pp collisions indepen- 
dently of the choice of actual jet algorithm or observable fl2l . 
While the formalism of SCET leads directly to expressing these 
observables as convolutions of separate hard, jet, and soft func- 
tions, blind use of this procedure without considering further 
specific properties of each chosen observable can hide whether 
factorization truly holds in a particular case or not. 
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where in the first form E, is the energy of particle i and f?, is the 
angle between its momentum and the thrust axis of X. In the 
second form, pf is the momentum of particle i transverse to the 
thrust axis, and 77, is its rapidity with respect to the direction of 
the thrust axis. We assume all final-state particles are massless. 

In a separate publication, using SCET, we calculate the an- 
gularity jet and soft functions to next-to-leading order in the 
strong coupling a s , resum large logarithms using renormaliza- 
tion group evolution, and model the nonperturbative soft func- 
tion in a way that avoids renormalon ambiguities IfTTIl . 

In this Letter, using angularity distributions as an example, 
we describe a simple, intuitive method for testing the validity of 
a factorization theorem deduced from the simple rules of SCET. 
We begin by naively presuming the factorizability of a given 
observable and then attempt to calculate perturbatively the one- 
loop jet and soft functions. If the factorization holds, each of 
these functions should be infrared-safe. If they are not, we learn 
immediately that the factorization breaks down. 

Perturbative infrared-safety of jet and soft functions is not, 
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of course, by itself sufficient to guarantee validity of the pro- 
posed factorization theorem. The size of power corrections 
must also be taken into account. The methods we describe in 
this Letter address only the former issue, not the latter. (Power 
corrections for angularity distributions and their implications 
for factorizability were studied in IfTUl [13] [18l .) However, our 
method is a quick and direct way to narrow down the class of 
observables for which a "generic" factorization deduced from 
SCET (e.g. D3) could actually be valid. 

Our analysis also sheds light on some issues related to in- 
frared divergences in effective theory loop integrals. Finding a 
tractable regulator in SCET that suitably controls all infrared di- 
vergences has been very challenging (see, e.g., lfT9ll2"0"ll ). Care 
is also required to define the effective theory such that it avoids 
double-counting momentum regions and infrared divergences 
of full theory diagrams. The procedure of zero-bin subtraction 
has been proposed to eliminate such double-counting lf2T1 . 

We will address each of these issues without explicit cal- 
culation of the jet and soft loop integrals or use of an explicit 
infrared regulator. Instead we just examine the regions of in- 
tegration surviving in the sum over all relevant diagrams. We 
will work in pure dimensional regularization, and learn how to 
identify 1/e poles as infrared or ultraviolet in origin, clarifying 
the contribution made by scaleless integrals which are formally 
zero. We will thus conclude that the analysis is independent 
of the choice of any explicit IR regulator. In the process, we 
demonstrate the crucial role of zero-bin subtractions in obtain- 
ing physically-consistent, infrared-safe jet functions in angular- 
ity distributions for all a < 1. The ideas and methods illustrated 
through our discussion of angularity distributions are more gen- 
erally applicable to other observables as well. 

2. Angularity Distributions in SCET 

The factorization theorem for the angularity distributions 
dcr/dT a takes the form, 



1 da 
cr dr, 



T - = H(Q; fx) j dr" a di» a dr s a 5{r a -T» a -T* a - t£) 
x/2(T>)7>>)S a (T>), 



(2) 



where cr is the total e + e~ — > qq Born cross-section, H is a 
hard function given in the effective theory by the square of a 
matching coefficient dependent only on short-distance effects, 
are jet functions dependent on the partonic branching and 
evolution of each of the two back-to-back final state jets, and S a 
is a soft function dependent on the low energy radiation from 
the jets and the color exchange between them. All the func- 
tions depend on the factorization scale p, with this dependence 
cancelling in the full cross-section. The factorization theorem 
Eq. ([2]) for angularity distributions has been proved in full QCD 
Ifl3ll and in SCET |fT0l[18l , for a < 1, where this condition was 
derived from the size of power corrections induced by replac- 
ing the thrust axis implicit in Eq. ([I} with the collinear jet axis n 
iTTOl n~3l . Our attempt to calculate perturbatively the jet and soft 
functions in Eq. (|2) will provide a complementary way to de- 
duce this condition and an intuitive explanation of the absence 



of infrared divergences in the jet and soft functions for a < 1 
and their appearance for a > 1 . 

Collinear and soft modes in SCET are distinguished by the 
scaling of the momenta of the particles they describe. The 
light-cone components p = (n ■ p, h ■ p, p x ) of collinear modes, 
where n,h are light-cone vectors in the +z directions, scale as 
Q(A 2 , I, A) or Q(l,A 2 , A), and soft modes as Q(A 2 , A 2 , A 2 ). Q is 
the hard energy scale in the process being considered (here, the 
center-of-mass energy in e + e~ collisions), and A is a small ratio 
of energy scales, here A = ^Aqcd/Q- Collinear momenta p c 
are split into a "label" piece p c containing the order Q and QA 
momenta, and a "residual" piece k c all of whose components 
are order QA 2 . A redefinition of the collinear fields through 
multiplication by soft Wilson lines decouples soft and collinear 
modes in the SCET Lagrangian to leading order in A fl6). 

The soft function S a in Eq. |2]) is defined by 

S a {T s a -p) = T^Tr<0| Yl(0)Yl(O)S(T s a - f s a )Y n (0)Y h (0) |0> , (3) 



and the jet functions by 



2 I ' a/3 N C 



L Tr f>*e a * (4a) 

X {0\Xn,Q(x)aS(T" a - OP».e(0)/J |0> 

y "«-" ) (*t=^ Tr /^/ rf4 -"'"' <4b) 

x (OlXn.-QWM - r a )Xn,- Q (0) a |0) . 

The traces are over colors, the light-cone momenta are defined 
l + - n ■ I and ]C — n ■ k, and the subscripts Q on the jet fields 
in Eq. Q specify that they create jets with total label momenta 
Qn/2 and Qn/2 \5\. The soft Wilson line Y„ in the soft function 
is the path-ordered exponential of soft gluons, 



T„(z) = ,Pexp 



f 

Jo 



dsn • A s (ns + z) 



(5) 



and similarly for Y„, with the bar denoting the anti-fundamental 
representation. The fields Xn,n m the jet function are the product 
of collinear Wilson lines and quarks,^-,, = W„£„, where 



W n (z) = Pexp 



I dsn ■ A„(h 
Jo 



s + z) 



(6) 



and similarly for Wn- The operator f a acts on final states \X) 
according to 



fa\X) 



1 £| pf | e -M(>-.)| X> 

ieX 



(7) 



and is constructed from the energy-momentum tensor T^ v [22 1, 
and the operators f"'"' s in Eqs. ([3]) and Q are constructed by 
keeping only the n, n-collinear or soft terms in T^ v [10|. For 
further details of the SCET Lagrangian and the Feynman rules, 
we refer the reader to Refs. lH|5]|6l. 

Next we proceed to examine the infrared behavior of 0(a s ) 
contributions to the soft and jet functions of Eqs. ([3| and Q 
calculated in perturbation theory. 
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3. Divergences in the Soft Function 

The soft function is calculated from the cut diagrams in 
Fig.[T] with an additional delta function 6(r s a - T a (X s )) inserted 
along the cut, where X s is the final state created by the cut, and 
T a (X) is given by Eq. ([1]). This modified cutting rule is required 
by the insertion of the 6{T s a - f operator in Eq. ([3) [ 17 1. 

In diagrams (A) and (B) of Fig. [T] with a virtual gluon, this 
delta function is just 6(j s a ), whose coefficient is given by the vir- 
tual gluon loop integral. Using pure dimensional regularization 
in d = 4 - 2e dimensions, this integral is scaleless and de- 
fined to be zero. This zero is actually a quantity proportional to 
1/euv - 1/fiR, and ordinarily plays the role of cancelling 1 /eiR 
divergences in diagrams (C) and (D) in which the cut creates a 
real gluon, and converting them to 1 / euv Il23l |2D . 

It is not at all obvious, however, how this cancellation can 
occur, since the virtual diagrams are independent of a while 
the real gluon diagrams depend explicitly on a. One often just 
prescribes the virtual diagrams to take a form that converts the 
l/e IR poles in the real diagrams to 1/euv, but this prescription 
is ad hoc and, as we will see below, potentially misleading. 

The soft function takes the general form 



S a (T s a )=A6(T s a ) + Y J B n 



0(T*)log"T* ( 



(8) 



Since the virtual diagrams are proportional to 6{T s a ), to study 
how they cancel the IR poles in the real diagrams, we only need 
to isolate the coefficient A of 6(rl). By integrating all the dia- 
grams over between and 1 , using the property of the plus 

functions J dx[8(x)(log" x)/x]+ = 0, we isolate A. We will 
denote as /J and /* respectively the virtual and real diagrams' 
contributions to A. 

The virtual diagrams' contribution to A is 



's 



a s C F {An^Y 
~r(l-e) 



dk" 



dk-(k + k-y l - £ , 



(9) 



whose integration region is the entire first quadrant of the ^ 
plane as shown in Fig.|2j\. 

In the real gluon diagrams, the cut creates a state X g with 
a single soft gluon, and the operator 5(T n a - f") acting on X g 
introduces the delta function 6(T s a -T a (X g )) into the integral over 
the gluon momentum k, where 



Ta(X„) 



~ Q \ k ~ 



il-S 



i* + i s 



for 
for 



r >k + 

r <k + . 



(10) 




Figure 1: The (A), (B) virtual and (C), (D) real gluon contributions to the soft 
function. The gluons all have momentum k. 
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Figure 2: Regions of integration in the k~ , k + plane for the coefficient of 6(r s Q ) 
in the a = soft function Sq(t s q ). (A) The region of integration for the virtual 
diagrams is the entire first quadrant. (B) For the real diagrams the region is 
S, which contains IR divergent regions. (C) These are converted in the sum of 
virtual and real graphs into the purely UV region S. 
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Figure 3: Region of integration in the k" 
the soft function S„(t£) for a = -1, a = 



, k + plane for the coefficient of S(T s a ) in 
0.5, and a = 1. The region S is formed 



by summing real and virtual diagram regions as in Fig. [2] For a < l,S always 
remains above the line k + k~ = Q 2 , which is the boundary of the a = 1 region 
and divides the infrared and ultraviolet regions of the soft loop integration. 



The real diagrams thus contribute 

a^iA^y rr dk+dk - (k+n ^ 

n T(l - e) Sis 

to A. 1^ depends explicitly on a through the integration region 
determined by the delta function <5(t£— T a (X g )), restricting gluon 
momenta k to the region S, given by fe* > and 



(12) 



(r)3(Jt + ) I_ 5 < Q for k~>k + 
(k + y-(k-y-i < Q for k-<k + . 

This region is plotted in Fig.|2b for a — 0. 

The virtual and real integrals I s ' contain exactly the same 
integrand, but with opposite relative signs and integrated over 
different regions of the k* plane. Thus, in the sum of the virtual 
and real integrals /* + /J, the integrals over the overlapping part 
of the regions cancel, leaving an integral over the region S, as 
illustrated in Fig. |2p. 

In Fig. [3] we plot the integration regions resulting from the 
sum of virtual and real diagrams for several other values of a. 
For a < 0, the resulting region of integration >S always satisfies 
fc* > Q, and is manifestly a purely UV region. Between < 
a < 1, A:* do approach zero on the boundaries, but the product 
k + k~ ~ always is greater than Q 2 , so divergences in the soft 
loop integral are still purely UV. For a = 1, the boundary of 
the region is the line of constant k + k~ = Q 2 . For a > 1, the 
boundary drops below this line, so, as — > oo, the product 
drops to k+kr ~ k 2 ± — > in the region along the boundary. But 
such momenta are in fact collinear. The sum of soft diagrams 
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still contains contributions from collinear modes. Continuing 
to explicitly evaluate /| + /J, we find this sum of integrals is 
convergent for a < 1 when e > 0, so the poles are 1/euv, but 
not for a > 1, in which case uncancelled IR divergences remain. 

The shapes of the regions in Figs. [2] and [3] also tell us how 
using an explicit IR regulator would affect our analysis, and in 
fact teaches us that the choice of regulator must be made with 
care. For example, we might choose an effective cutoff A on fc* 
in soft loop integrals as used in E4l . in which the soft function 
for jet energy distributions was calculated to one loop and ar- 
gued to be IR finite. In these cases this regulator successfully 
cuts off the divergences arising from the regions ^ — > 0. How- 
ever, using this regulator for < a < 1, we find that the soft 
function still contains log A divergences and 1 / euv divergences 
even though the above analysis shows that + /J is actually 
IR finite. From Fig. [5]it is evident that a lower cutoff on ]& also 
cuts off regions where fe* — > oo, so it acts also partially as a UV 
cutoff. This underscores the challenge of defining consistent, 
explicit IR regulators in SCET lfT9ll20l . 

We draw two lessons from the analysis thus far. The first is 
that in pure dimensional regularization, the coefficient of (1/euv - 
1 / 6ir) in a virtual diagram cannot be determined from the vir- 
tual diagram alone, but only together with the real diagram 
whose IR divergence it is supposed to cancel (cf. G5ll ). The 
reason that the virtual subtraction can depend on a even though 
by itself it is independent of a is that the area of overlap between 
the integration regions of real and virtual diagrams depends on 
a. The second is that the presence or absence of IR divergences 
in the sum of all contributing loop integrals can (and should) 
be determined before completely evaluating the integral with a 
given IR regulator. Looking at the shape of the region of in- 
tegration in momentum space as above avoids confusion about 
the consistency of the regulator itself. 

4. Divergences in the Jet Function 

Now we analyze the jet functions Eq. Q at 0{a s ) in pertur- 
bation theory. We will observe the same breakdown of infrared 
safety as a — > 1 due to the momentum regions beginning to 
include IR divergent regions. We will consider just the jet func- 
tion J"(T" a ); identical analysis applies to J„(t"). 

The diagrams contributing to ./"(t") in Eq. ( [4a] ) are shown 
in Fig. [4] In graphs (A) and (B), the gluon is emitted from 
a collinear Wilson line W„ or W,J in the jet fields Xn>Xn- The 
sum of graphs (C) and (D) is equivalent to graphs in full QCD 
by a field redefinition |26 27 1 and is manifestly IR finite for 
all a < 2, and we will not consider them further here [17|. 
The total momentum flowing through each diagram is Qn/2 + 1, 



(A) (B) (C) 

Figure 4: Diagrams contributing to the angularity jet function J%(t%). The total 
momentum through each graph is Qn/2 + /, and each gluon momentum is q. 
(A) Wilson line emission diagram and (B) its mirror; (C) sunset and (D) tadpole 
QCD-like diagrams. 



with the label component Qn/2 specified by the labels on the 
jet fields in the matrix elements in Eq. ( |4"a| >, and I the residual 
momentum. The total momentum of the gluon in each loop is 
q — q + q r , where q is the label momentum and q r the residual 
momentum. The diagrams must be cut in all possible places, 
and a delta function 6{T n a - T a (X n )) inserted along each cut. To 
obtain J"(t^) we then integrate over l + according to Eq. Q. 

An integral over the collinear gluon momentum q is a sum 
over q and an integral over q r . The sum excludes the value 
q = 0. The sum and integral can be combined into a single con- 
tinuous integral over the total q if a "zero-bin subtraction" is 
also taken ET1 . which avoids the double-counting of soft con- 
tributions between the soft and jet functions lfT3l[T8ll20ll25ll28l . 
Below we will refer to integrals or graphs before the zero-bin 
subtraction as "naive", and after the subtraction "collinear". 

Virtual graphs and zero-bin subtractions in pure dimensional 
regularization again contain scaleless integrals, which are zero, 
but as we observed in the calculation of the soft function, we 
must examine the regions of integration to observe how the can- 
cellation of IR divergences among all the graphs occurs. 

Diagrams created by cuts through the single quark propa- 
gator in graphs (A) and (B) in Fig.|4]leave a virtual gluon loop 
and are proportional to S(t"), whose coefficient we extract. The 
naive virtual graph contributes 



~y _ a s C F {Ann 1 ) 



-e)Jo Jo Nqi) 1+e W- Q) 



(13) 

which goes over region 'VinFig.m The zero-bin subtraction 
from the virtual graph is 



r = - 



a s C F (4^-yt/ 2 ) 



r(i 



% fV f 

- e) Jo Jo 



dql 



(q 2 J 1+e q- 



(14) 



which goes over the whole first quadrant, region *Vo in Fig. [5] 
So the total virtual collinear contribution iY = IY - I v ,, is 



a s C F (4nfi 2 y 
n F(l - e) 



Jo q± (qi) 1+6 1 J Q q~ + Jo 



dq 

~Q 



(15) 

where the l/q~ term is integrated over region *V in Fig. [5] and 
the 1/Q term over *V. 

Now we add the contribution of the graphs (A) and (B) in 
Fig.|4]cutting through the gluon loop, creating a final state X qg 
with a collinear quark and gluon, with 



T a(X q „) 



^[(9 _ ) ! (9 + ) 1_! 



Q 



+ (Q-q-)nr 



| (16) 



and insert £(r" - T a (X qg )) into the integral over the gluon mo- 
mentum q. As in the case of the soft function, we need only to 
isolate the coefficient of (5(t") in the jet function J%(t%) to study 
the cancellation of IR divergences with the virtual graphs. We 
do so by again integrating over < t„ < 1 . 

The contribution to the coefficient of <5(t^) from the naive 
Wilson line graphs (A) and (B) in Fig.|4]with a cut through the 
gluon loop is 



JR 



a s C F (4JT/J 2 ) 



n T(l 



dq dq]_ 
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(17) 
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Figure 5: Regions of integration in the q~,q 2 ± plane for virtual gluon diagram 
contributions to the coefficient of 5(r") in the jet function J'^(r"). 'V is the 
region for the naive integral, *Vo for the zero-bin subtraction, and 'V for the 
sum of these two contributions. 




Q Q Q q 

Figure 6: Regions of integration in the g~,q^ plane for real gluon diagram 
contributions to the coefficient of <5(rjJ) in the a = jet function 7Jj(rJJ). ft is the 
region for the naive integral, ftg for the zero-bin subtraction, and ft for the sum 
of these two contributions. 

where % is the region in the first quadrant of the q~,q]_ plane 
under the curve 
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shown for a = in Fig. [6] The zero-bin subtraction is 

R a s c F (M^y rr 2 1 1 

where Kq is the region given by q~ > and 
0<ql < [Q(q-y- a ]^ 2 , 



(18) 



(19) 



shown for a = in Fig. [6] Subtracting the two integrals, /* 
cti 

R _ a s C F (47Tfi 2 y 



(20) 

R _ 



- /* Q , yields the correct collinear integral 



n F(l-e) 



1 1 



(21) 



where the first integral goes over the region % formed by re- 
moving % from 7?o, illustrated in Fig. a — 0. The second 
integral, containing l/Q, still goes over %. 

Thus, the sum /„ = /J + of the collinear virtual and real 
graphs is 



/„ = - 



a s C F (47in 2 y 



n F(l - e) 
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where the regions are shown in Fig.|7]for several values of 
a, and the notation JTi - J2 means the integrand has an extra 



91 
4 



a = -2 









a = 









4 



a = 0.9 


J2 


Ji 




— — , 





a = 1 




Ji 









Q 



Q 



Figure 7: Regions of integration formed by combining real and virtual diagram 
regions in Figs.|5]and[6] The regions J~i 2 result from subtracting the virtual 
diagram region *V in Fig.[5]from the real diagram region ft in Fig. [6] For a = 1, 
the region encounters an unregulated divergence at q~ = 0. 



minus sign in contains only UV divergences. For all 

a < 1, STi avoids the boundary at q~ = 0, and the integral is 
convergent for e > 0. The 1 /e poles in this integral, as well as 
in the integral on the second line of Eq. p2t, are then purely 
UV. For a — 1 , however, the region Ji reaches the boundary 
at q~ = 0, and the integral is no longer finite. The jet function 
is not infrared safe for a > 1, just as we found for the soft 
function. 

Although the full distribution dcr/dT n is infrared safe for 
a < 2, for a > 1, contributions of the soft and collinear modes 
of SCET with the momentum scalings specified in Sec. [2] (so- 
called SCETi modes) do not entirely separate from each other. 
The soft integration regions illustrated in Fig.[3]for a < 1 grow 
for a > 1 to include the contribution of collinear modes, and the 
collinear integration regions in Fig. [7] grow to include modes 
which are soft. Angularity distributions with a > 1 are domi- 
nated by jets so narrow that collinear and soft modes have the 
same virtuality of order Aqcd- We observe this in a full calcula- 
tion of the jet and soft functions to 0(a s ) ifTTl . which manifests 
the natural scales in the jet and soft functions where large log- 
arithms are minimized, /uj = QtJ a ^ and fi$ = Qr a , which 
become equal at a = 1 . Thus, the separation of scales required 
by the formalism of SCETi no longer holds. In this case, the 
modes may be distinguished by their rapidity, as was proposed 
in the formalism of SCET n ED . 

5. Conclusions 

Although hard-scattering cross-sections can be written for- 
mally in a factorized form based on naive application of a for- 
malism such as SCET, the properties of the chosen observable 
determine whether or not the effective theory is applicable and, 
so, whether the factorization theorem is actually valid. Such 
a theorem must pass a number of tests. The method we have 
presented is a straightforward and intuitive test of the infrared- 
safety of jet and soft functions in a proposed factorization theo- 
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rem. We illustrated the method at 0(a s ) for angularities, whose 
tunable parameter a allowed us to study the continuous progres- 
sion from infrared-safety of jet and soft functions for a < 1 
to its breakdown for a > 1, but the method is more gener- 
ally applicable to other observables as well. The test will re- 
veal those observables for which the naive SCETj factorization 
fails. Through our analysis, we have illustrated the crucial role 
of zero-bin subtractions in effective field theory, and the man- 
ner in which scaleless integrals in pure dimensional regular- 
ization convert IR into UV divergences in infrared-safe quanti- 
ties, without choosing any ad hoc prescriptions, allowing one 
to classify IR and UV divergences independently of an explicit 
IR regulator. 
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